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Topological phases of matter that depend for their existence on interactions are fundamentally interesting and 
potentially useful as platforms for future quantum computers. Despite the multitude of theoretical proposals the 
only interaction-enabled topological phase experimentally observed is the fractional quantum Hall liquid. To 
help identify other systems that can give rise to such phases we present in this work a detailed study of the effect 
of interactions on Majorana zero modes bound to vortices in a superconducting surface of a 3D topological 
insulator. This system is of interest because, as was recently pointed out, it can be tuned into the regime of 
strong interactions. We start with a OD system suggesting an experimental realization of the interaction-induced 
Zs ground state periodicity previously discussed by Fidkowski and Kitaev. We argue that the periodicity is 
experimentally observable using a tunnel probe. We then focus on interaction-enabled crystalline topological 
phases that can be built with the Majoranas in a vortex lattice in higher dimensions. In 1D we identify an 
interesting exactly solvable model which is related to a previously discussed one that exhibits an interaction- 
enabled topological phase. We study these models using analytical techniques, exact numerical diagonalization 
(ED) and density matrix renormalization group (DMRG). Our results confirm the existence of the interaction- 
enabled topological phase and clarify the nature of the quantum phase transition that leads to it. We finish with 


a discussion of models in dimensions 2 and 3 that produce similar interaction-enabled topological phases. 


I. INTRODUCTION 


Non-interacting electron systems have recently provided a 
new playground for experimental and theoretical research — 
topological insulators (TIs) and topological superconductors 
(TSCs) Hie The hallmarks of these materials are symmetry- 
protected gapless edge modes. Both fundamental and crys- 
talline symmetries can play a role in protecting these states “© 
The TIs and TSCs are extensively studied both on their 
own and as building blocks in engineered exotic quantum 
states /714) 

Majorana bound state is an example of such an edge state, 
which can emerge both at the edge of a 1D TSC and in 
the vortex of 2D TSC, which may be either intrinsic or re- 
alized at the interface between a 3D TI and an ordinary 
superconductor 205] Such Majorana bound states occur at 
zero energy and are thus often called Majorana zero modes 
(MZMs). The effort to detect such MZMs has been focused 
on the zero-bias conductance peak /iH8) indicative of the den- 
sity of states at zero energy, which is, however, not a defini- 
tive proof of a Majorana bound state! The physics of the 
MZMs still to be observed is much richer: most importantly, 
upon adiabatic exchange they behave as non-abelian anyons 
and provide a pathway to (non-universal) topological quan- 
tum computation 22 

It is known that the electron-electron interactions enrich the 
palette of topological states. The first discovered topologi- 
cal material with inherently strong electron-electron interac- 
tions is the fractional quantum Hall effect 2223] For certain fill- 
ing fractions the emergent quasiparticles are predicted to be 
Majoranas, or even more exotic Fibonacci anyons, that may 
provide a platform for universal quantum computation!22# 
These systems require extreme conditions to be probed and 
used: very clean samples, low temperature, and high mag- 
netic fields. One way to avoid these complications is to en- 
gineer the Majorana bound states2524 and use the interactions 


between them to construct universal quantum computer. Ma- 
jorana surface codes realizing exotic non-abelian quasiparti- 
cles allowing for universal quantum computation, have been 
theoretically proposed 


For the experimental realization of these proposals there are 
several preliminary steps to be made. Though there have been 
indications of the Majorana bound states in the conductance 
measurements in several systems, the possibility to induce 
interactions between them still remains unconfirmed. In the 
present paper we endeavor to fill this gap from the theory side 
by suggesting realistic setups that can be used to probe the 
interactions between Majoranas. 


As a basis we use a system proposed recently as a conve- 
nient platform to study strong interactions between Majorana 
zero modes: vortices in a surface of 3D TI with induced su- 
perconductivity and the chemical potential tuned to the Dirac 
point. It was argued that for this special point direct Majorana- 
Majorana tunneling is absent due to the additional chiral sym- 
metry of the underlying physical system 28131] The system of 
MZMs is then in the so-called symmetry class BDI and can 
be viewed as respecting a fictitious time-reversal symmetry=2! 
© such that ©? = 1. The dominant term in the low-energy 
Hamiltonian is the 4-Majorana interaction that arises from 
the Coulomb interaction between the constituent electrons!) 


A number of interesting phases have been identified in this 
setup2#34 including some interaction-enabled phases 22132139] 


In the present paper we study the effects of interactions on 
the phases of the MZMs in dimensions 0-3. We consider sys- 
tems with time-reversal symmetry ©, and the main focus of 
the paper is on OD and 1D. We start with the OD case. We 
first argue for Zg Fidkowski-Kitaev periodicity of the ground 
state degenerac to coincide with the degeneracy of the 
entanglement spectrum 23 The problem can be stated as fol- 
lows. Given the set of n MZMs such that no direct tunneling 
between them is allowed, what is the ground state degener- 
acy in the presence of generic 4-fermion interactions? By a 


combination of simple arguments and direct computations we 
determine the degeneracy for n = 1...8 and then present an 
argument that the pattern repeats itself with the periodicity of 
8. We then suggest a tunneling probe as a way to experimen- 
tally observe the expected Zg periodicity. The probe would 
reveal the absence of the zero-bias peak for the total vorticity 
being a multiple of 4, in contrast with the non-interacting case 
for which we expect the zero bias peak to be absent for even 
total vorticity. 

The second part of the paper focuses on the interaction- 
enabled crystalline topological phases in dimensions 1, 2 and 
3. As previously noted by Lapa, Teo and Hughes (LTH) 4! 
in the presence of inversion and the fictitcious time-reversal 
symmetry there exist no topological phases in non-interacting 
1D models. This is because the inversion symmetry maps 
the integer topological invariant v to —y. For an inversion- 
symmetric system, therefore v = —v, implying vy = 0 as 
the only solution. In the presence of interactions the integer 
classification changes to Zg. Under Zg, remarkably, the equa- 
tion v = —v has a non-trivial solution v = 4, indicating a 
genuine interaction-enabled crystalline topological phase. We 
discuss a physical setup that realizes such a phase in the sys- 
tem of vortices and antivortices in the surface of a 3D TI. 
We start this discussion by introducing a non-trivial but ex- 
actly solvable “2-leg LTH model” that is closely related to 
the LTH model but does not exhibit the topological phase. It 
has two distinct phases, distinguished by a conventional bro- 
ken symmetry. The 2-leg LTH model can be mapped onto a 
collection of the Kitaev chain models, which allows for ex- 
act solution. We use this model to benchmark our numerical 
simulations which we then employ to study the “4-leg LTH 
model” which has an interaction-enabled topological phase. 
The 4-leg model is exactly solvable only in two extreme lim- 
its, the strong coupling limit shows the topological phase. We 
use ED and DMRG techniques to study its ground state de- 
generacy, central charge and the entanglement spectrum away 
from the exactly solvable limits. We show that the topological 
phase extends over a significant part of the phase diagram. 

We close by proposing variants of the interaction-enabled 
topological crystalline phases in 2D and 3D that can be built, 
at least in principle, from the ingredients at hand. As in 1D, 
these phases require interactions for their existence and ex- 
hibit anomalous edge or surface states protected by a combi- 
nation of time-reversal and crystalline symmetry such as in- 
version or discrete rotation. 


Il. PHYSICAL SYSTEM 


In this section we introduce our platform for strongly inter- 
acting Majorana fermions and review some facts that will be 
important in what follows. We consider vortices in the sur- 
face of a strong topological insulator in the proximity with 
an s-wave superconductor, which was first proposed by Fu 
and Kane!5| When the chemical potential is tuned to the Dirac 
point a fictitious time reversal symmetry © emerges. Hop- 
ping of MZMs located in the individual vortices is forbid- 
den by © and the second major effect, interaction, becomes 
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dominant®22/ Here we will only review the surface physics 
briefly and point out the important for the present work prop- 
erties of the Hamiltonian and the wavefunctions in the Fu- 
Kane model. The Hamiltonian of such a system in the Nambu 
four-vector basis (Ctr, Cyr, a —cl.) reads 


Hypx(r, P) =UP2OnTz + VPyTyTz — PT x 
+ ReA(r)tz + ImA(r)t,, (1) 


where o and 7 are Pauli matrices in the spin and Nambu ba- 
sis respectively. Particle-hole symmetry automatically is pre- 
served with = = t¥a0¥%Kk, where K represents complex con- 
jugation. In the absence of the vortices the system also pre- 
serves physical time-reversal symmetry, 0 = o, kK. By fine 
tuning of the system, namely adjusting jz to zero, the system 
respects time-reversal symmetry with 9 = t*a*K. Notice 
that the physical time-reversal symmetry, which is broken in 
the presence of magnetic field, squares to —1, as it should 
for spin-4 particles, but the ©? = 1, which is allowed for an 
accidental symmetry. Furthermore, chiral symmetry is auto- 


matically preserved with the operator I] = =O = —r’o”%. In 
the following we focus on the j = 0 situation with the extra 
symmetry. 


If we assume a vortex at the origin, the gap profile in polar 
coordinates is A(r) = A(r)e’”®, where n is the vorticity of 
the vortex. For a system respecting © and E the vorticity is the 
same as the total number of the MZMs inside the vortex as we 
will see shortly. We call the cases with n > 0 vortices, and 
n < 0 anti-vortices. The wavefunctions for these two cases 
are respectively 


B, =(0,u,uj,0)7, B= (uy,0,0,-uk)” (2) 


These wavefunctions belong to the opposite eigenvalues +1 of 
the chirality operator II. Importantly, the symmetry II forbids 
the hybridization between two Majorana zero modes with the 
same chirality28!31) 

To simplify the Hamiltonian and bring it into block off- 
diagonal form we perform a unitary transformation Hp — 
UHpKU~! with 


0 
1 
o|- (3) 
0 


This renders the equations for the zero-mode for both vortici- 
ties 


D(r)x (tr) =0, Di (r)xy = 0, (4) 
where 
e irp-ia Ay (r) 


ae (-i9, - 20 ) 
ele (—i9, + 22) 


eine ria No (r) 


i= (5) 


and &; = U(0,x,)7, ®+ = U(x+4,0)". We thus establish 
the basic building block of the setups we will consider — a 
vortex of vorticity n, which carries n Majorana zero modes. 


If n = 1, then the zero mode obtained from the equation (4) is 
automatically the Majorana zero mode, otherwise MZMs can 
be obtained as linear combinations of these modes. 

In what follows we shall call a; the creation operators of 
the MZMs from the | sector and 8; from the + sector. We 
see from the explicit form of the Majorana wavefunctions 
in the two spin sectors, see Eq. 2). that the two types of 
MZMs transform differently under the emergent time rever- 
sal, namely=2! 

0a,;O-! => as, 05,071 = —B;. (6) 
Along with the property 010~-! = —i this implies that one 
can construct a complex fermion c; = $(a; +i8,) that trans- 
forms naturally20#144) under @, as Oc; e-l= c;. In the rest 
of this paper we shall study systems of MZMs described by 
Q-invariant Hamiltonians that contain both 2-fermion “tun- 
neling” terms and 4-fermion interaction terms. The former 
arise from the wavefunction overlaps between MZMs while 
the latter originate from the interactions between the con- 
stituent electron degrees of freedom. These include the direct 
Coulomb repulsion, as well as interactions mediated by other 
degrees of freedom such as phonons, which may be attractive. 

In a O-invariant Hamiltonian Eq. (6) allows tunneling terms 
between different types of MZMs of the form 


Hin = is? tijaiB;, (7) 


but prohibits those between the same chirality such as ia;a; 
and 76;,;. Interaction terms with an even number of a oper- 
ators such as Qj;ajAKAm, BiB; Bem and aja; b~Bm are al- 
lowed but those with an odd number are prohibited. An in- 
teresting consequence of this structure arises in a system that 
is composed solely of vortices (and no antivortices). Accord- 
ing to the above discussion all tunneling terms are then for- 
bidden and the kinetic energy is quenched. Interaction terms 
however are allowed by symmetry and the system is, there- 
fore, inherently strongly interacting. In a realistic setting, the 
chemical potential will always be slightly detuned from zero 
(either globally or locally due to disorder) leading to non-zero 
Hin. However, as long as this detuning is small compared to 
the interaction scale g the system must be viewed as strongly 
interacting. 

In Sec. II below we work with interactions only (assuming 
jt = 0 and the presence of only vortices in the system) while 
in Sec. IV we consider interacting systems composed of both 
vortices and antivortices that exhibit interactions as well as 
hopping. In this way we establish experimentally realizable 
strongly interacting phases of the Majorana bound states. 


Il. FIDKOWSKI-KITAEV PERIODICITY 


In their seminal work Fidkowski and Kitae#!l predicted 
that in the presence of interactions the Z topological classifi- 
cation of a 1D systems of fermions with symmetry © (class 
BDI) breaks down to a periodic Zg one. That is, the gapped 
phase with the winding number v = 8 can be continuously 
deformed to the trivial phase by turning on local interactions 
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FIG. 1: Top view on the setups we suggest to observe the Zg pe- 
riodicity of the BDI Majorana model. Blue is the area covered by 
the superconductor, and yellow is the bare surface of the topolog- 
ical insulator. a) A single large hole in the superconductor. Simi- 
lar to Corbino disk geometry all the magnetic field goes through the 
hole. Total flux is n®p and the total number of MZMs in the non- 
interacting model is n. b) The lattice of vortices. We consider a 
square lattice with 2 x n/2 vortices. If n is odd, there are (n — 1)/2 
lines of 2 Majoranas and one line with just one, similar to the n = 5 
case in the figure. Such an arrangement can be created by introducing 
strong pinning centers into the superconductor. 


without any gap closing. The original non-trivial phase pos- 
sesses 1 Majorana zero modes at each edge. The above result 
fundamentally depends on the fact that symmetry-preserving 
interactions can completely remove the 2*-fold ground state 
degeneracy of 8 MZMs but cannot do it for any smaller num- 
ber of MZMs. Fidkowski and Kitaev considered a very spe- 
cific, highly symmetric, form of interactions to demonstrate 
the above effect. How general or finely tuned the interactions 
need to be to produce the unique ground state for 8 MZMs 
however remained unclear. Here we show that under quite 
general conditions the ground state degeneracy of 8 Majo- 
ranas is split by Coulomb interactions between the constituent 
electrons. Moreover, we give a general argument showing the 
Zs periodicity of the ground state degeneracy in a system of 
n Majoranas in accordance with the Fidkowski-Kitaev result. 
We also explain how the resulting Zg periodicity induced by 
interactions can be experimentally observed. 


In this section we are interested in compact OD structures 
with the total vorticity n. We assume that 4 = 0 and no tun- 
neling terms between MZMs are therefore allowed although 
generic 4-fermion interaction terms are present. We consider 
two specific physical setups realizing such a system, shown 
in Fig. |1| The first is a superconductor with a large hole put 
on top of the TI with n superconducting flux quanta threaded 
through the hole as outlined in Fi g.[Ifa). This geometry can be 
viewed as a Corbino disk familiar from many superconduct- 
ing applications. We therefore expect that this configuration 
is the easiest for experimental access. The second setup is an 
array of single vortices illustrated in Fig. [I{b). This requires 
more fine-tuning, for example, creating strong pinning centers 
for the vortices, but is ultimately also accessible experimen- 
tally, especially since our qualitative results do not depend on 
the exact geometric arrangement of vortices. We discuss the 
tunneling conductance signature of the interacting system and 
how it is modified as compared to the non-interacting case. 


A detailed calculation of the energy spectra in these two 
systems shows that the same ground state periodicity of 8 in 
n is observed in both. We generalize this result by providing 
an argument which proves the connection between the peri- 
odicity of the entanglement spectrum and the periodicity of 
the ground state degeneracy for a system of n Majoranas with 
generic 4-fermion interactions. We also comment on the ap- 
plicability of this argument to the physical realizations of the 
interacting Majorana model suggested above. 


A. Giant vortex in a Corbino geometry 


We start our discussion by considering the setup depicted 
in Fig. [Ifa). We envision a hole in the superconducting coat- 
ing, through which a number n of superconducting flux quanta 
®o = hc/2e is threaded. This is equivalent to systems already 
studied in the context of TIs#2!In this Section we will consider 
one sign of vorticity, positive for definiteness. We search for 
the solution of Eq. (4) in the form 


1 etme nls) ge (p) 
Xm(r) = Va ( e~impta/2—n/4)y, (7) » (8) 


to obtain 


Ao(r)tm(r) + (0; — @) um(r) = 0, 
{ Ao(r)¥m(r) + (Op — 2==2) um(r) =0. 


Let us now take into account the practical realization of this 
giant vortex. We imagine a hole drilled in the superconductor, 
thus leaving the order parameter zero inside the hole and large 
outside, see Fig. [I{a). Thus inside the hole the Eq. (9) reduces 
to 


(0, - nm) Um(r) = 0, 
{ (0, = ud) Um (Tr) = 0, (10) 
having the obvious solutions 
Um(r) x r™, (11) 
Wmir) o pt in™, (12) 


Notice that when either of the powers is negative, the solution 
becomes non-normalizable at r = 0, similar to the arbitrary 
gap profile situation discussed first by Jackiw and Rossi/“® 

The boundary condition at the external radius ro of the hole 
is determined by the absence of solutions growing into the su- 
perconductor. The equation (9) under the assumption of large 
gap inside the superconductor gets simplified to 


{ Ao(r)um(r) + OrUm(r) = 0, (13) 


This requires 


Um(r0) = Um(T0). (14) 
Therefore inside the hole we have 
Um(r) = A(r/ro)™, (15) 


A(r/ro)”1-™. (16) 


Um(r) = 


Normalization of the wavefunction then requires 


To 
[rare +e) =1, a7 
0 
which finally gives 
1 1 —1/2 
A=71,? 1 
"o Ee In —2m—1 ee 


We have thus found the wavefunctions of the Andreev 
bound states in the hole at zero energy. To proceed in a uni- 
form fashion we change to the Majoranas basis. For that we 
first rewrite the solutions of the BdG equation in the second 
quantization form 


1 ; 
Xm =e f de (oO metoMunlrely 
+e7tneta/2y,, (roe) (19) 


We notice that the connection between Andreev bound state 
with indices m and n—1—™: ¥m = Xn_1~m A Un—-1—-m = 
Um. Thus form < (n—1)/2 we can define MZM operator as 


Xm + Xn-1—m 


Am = 20 
5} (20) 

and form > (n — 1)/2 
Qn = jXm — Xn-1-—m (21) 


V2 
For m = (n — 1)/2 and n odd we also have 
Am = Xm- (22) 


The interaction between MZMs arises from the underlying 
electron-electron interaction 


1 
Hine = 5 / drdr'p(r)V(r,r)o(r’), (23) 
where p(r) = >>, ch,cor represents the electron density op- 


erator. The interaction matrix element between MZMs can be 
computed by expressing the charge density operators in terms 
of the eigenstates of Hp and projecting onto the zero-energy 
Majorana subspace, as described in Ref.(31] One then obtains 
MZM interaction Hamiltonian of the form 


Hint = S JijklAjiajakal. (24) 
ijkl 


The matrix element g;;,; is computed numerically from the 
eigenfunctions Um, Um and the electron-electron interaction 
potential V. In what follows we will use screened Coulomb 
potential 


Vo 
|r —r’| 


Vir,r’) = eT FoF l/rs (25) 


where Vo is the Coulomb interaction strength, and r, is the 
screening length. We use r; — oo below, but we have 


n mod 8 [O/1/2|/3|4/5|6|7 
GS degeneracy | /1| /2| V4| V8] V4| V8/V/4| v2 


TABLE I: Ground state degeneracy in the presence of generic 4- 
fermion interactions as a function of the total number of MZMs n 
modulo 8. Irrational ground state degeneracy for odd numbers of 
MZMs implies that if we bring an additional MZM into the system 
but do not couple it to the system of interest, then the ground state 
degeneracy is multiplied by \/2 and becomes rational. This is always 
a valid interpretation since in any physical systems the total number 
of MZMs is even. 


checked that the finite screening length does not change the 
results qualitatively. We will work in the units where the ra- 
dius of the hole is 1. In these units if the hole has a radius 
of 100 nm and taking into account the dielectric constants of 
the TI materials, ranging from e = 20 for HgTe to ~ 200 for 
BiyTe3 J interaction strength is Vo = 0.72 — 0.072 meV, and 
the magnetic field to create vorticity 8 in the core is © 65 mT. 


By exact numerical diagonalization of the manybody 
Hamiltonian Hj; with all 4-fermion interactions as computed 
above we find the energy spectra of the system for n up to 16. 
The low-lying eigenenergies for n up to 8 are displayed in 
Fig. 2fa). The structure here is easy to understand for small 
n. For n = 1,2,3 no 4-fermion term can be constructed, 
Hint = 0, and all the states are at zero energy. For n = 4 
we have Hint = gQ1Q203Q4 which results in a 2-fold de- 
generate ground state and 2-fold degenerate excited state. For 
n > 4 progressively larger number of interaction terms can be 
constructed and the energy levels must be found numerically. 
From this we deduce the ground state degeneracies and find 
that they repeat in n with periodicity of 8 as predicted by Fid- 
kowski and Kitae41 The degeneracies are listed in Table 
[I|and agree with degeneracies of the entanglement spectrum 
for the corresponding topological phases predicted in Ref.|43] 
This is true for a range of parameters, in particular we verified 
this result for different screening lengths r,. Based on this 
observation we hypothesize that the periodicity of the ground 
state degeneracies holds for larger vorticities and for a generic 
arrangement of vortices. We will prove this statement below 
and specify more carefully the conditions under which it is 
valid. 


The Zg periodicity is experimentally observable in a sin- 
gle electron tunneling spectroscopy. Similar to the non- 
interacting case we consider tunneling into the states of the 
nm Majoranas, either using scanning tunneling microscopy 
(STM) or a normal contact attached to the hole region. The 
tunneling can only take place if the matrix element of «' and 
w between one of the ground states and the state into which 
we are trying to tunnel is non-zero. Otherwise one needs two 
particles to tunnel simultaneously, the process smaller in the 
tunneling amplitude, which we will neglect. We depict the po- 
sitions of the lowest excited levels along with the possibility 
to observe them in tunneling conductance in Fig. fa). 
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FIG. 2: Positions of energy levels and conductance peaks as a func- 
tion of the total vorticity n for the two setups depicted in Fig. [I](a) 
and (b) correspondingly. Red cross depicts the position of the energy 
level observable as Lorenzian peaks in tunneling conductance. Blue 
bar shows the energy level different from the ground state by an even 
number of fermions. Such level is not observable by single-electron 
tunneling. 


B. Array of vortices 


Another possible setup to observe the Zg periodicity is de- 
picted in the Fig. [I{b). We consider a square array of single 
vortices of size 2 x n/2. For odd n the last vortex is alone in 
the last line. In such a setup we need to make an assumption 
about the screening of the electron-electron interactions. As 
an example, we assume that all the vortices are closely situ- 
ated, and the amplitude of the interaction between any 4 Majo- 
ranas is the same, V,. This is an unnecessary assumption — we 
have obtained qualitatively same results for different screen- 
ing lengths — but it removes multiple interaction strengths and 
simplifies the expressions. 


The interactions, however, have gependence on the phase 
difference between the vortex cores! This is the dependence 
we take into account to obtain the positions of energy levels 
of the system. These, and their respective observabilities by 
single-electron tunneling are depicted in Fig. [2{b). We assume 
the possibility to tunnel into any of the vortices of the array, 
i.e. the tunneling probe much larger than the size of the array. 
Again, nothing qualitatively changes if we tunnel into a single 
generic vortex. 


C. General argument 


In this section we will discuss the ground state degeneracy 
for the system of n interacting Majoranas. Due to locality 
of the edge states, this system is equivalent to an edge of a 
BDI wire with the quantum number n. We will use the fol- 
lowing definition of the ground state degeneracy for a system 
of n MZMs: it is the degeneracy of the ground state in the 
entire multi-dimensional parameter space of the Hamiltonian 
with all the interactions allowed by symmetry turned on. The 
degeneracy can be higher on a measure-zero regions of pa- 
rameter space, where the accidental degeneracies are present 
on top the obligatory ones. The simplest example of such a re- 
gion is the point where all the interactions are turned off. We 
assume there are no two regions of parameter space that are 
having non-zero measure and different ground state degenera- 
cies. This assumption is justified by the zero-dimensional na- 
ture of the system, where there is no locality condition which 
can prevent two states from splitting. 

It turns out that it is convenient to prove a more general 
statement: generically all the excited states in such systems 
will have the same degeneracy as the ground state. We will 
use the results of the more general argument and apply them 
to the ground state degeneracy. It is straightforward to expand 
the argument to all the excited states as well. 

We can explicitly construct the states with the ground state 
degeneracies as shown in Table [I|for k =0...7MZMs. Ex- 
amples of such constructions were just shown in the previous 
section. In all such examples the degeneracy of the ground 
and excited states was found the same as expected accord- 
ing to the general periodicity of Table [I] We now add to the 
group of & Majoranas m blocks, each containing 8 Majoranas 
with non-degenerate ground states. We thus obtain a system 
of nm = 8m+k Majoranas with the ground state degener- 
acy the same as for & interacting Majoranas. Coupling these 
blocks together with a small enough coupling cannot increase 
the degeneracy of the ground state, since one needs finite per- 
turbation to close the gap. This limits the ground state degen- 
eracy for n Majorana bound states from above to the values in 
Table[]] 

What is left is to prove that the above limit is actually the 
ground state degeneracy for n Majoranas in a region of pa- 
rameter space. For k = 0,1,7 the limitation of the ground 
state degeneracy from below is obvious: for k = 0 the degen- 
eracy is 1, while for k = 1 and 7 it is \/2. The ground state 
degeneracy of at least \/2 is always present for an odd number 
of Majoranas since then there is a single Majorana mode that 
can be decoupled from all the others, thus giving the desired 
\/2 degeneracy. 

We will prove the limitation on the ground state degener- 
acy from below by induction. The starting point is that for 
n = 0,1, 2,3 no interactions term can be added to the Hamil- 
tonian, thus the limitation of the ground state degeneracy from 
below is obvious. For n = 4 there is a single interaction term 
allowed, as discussed above, it involves all the 4 Majoranas. 
Such an interaction halves the ground state degeneracy, and 
makes a doubly-degenerate excited staté22| We now come to 
the step of induction, which will generalize the result to arbi- 


trary n. 

The induction step goes via reductio ad absurdum path. As- 
sume that we have for n = 8m + k Majoranas degeneracy 
lower than expected. 

We already know that this is impossible for k = 0,1, 7. 
Let us proceed with the k = 2,6 case. The ground state de- 
generacy expected one reads off Table[l]is 2. Our assumption 
dictates that there is a region of the parameter space where 
the degeneracy is 1 (and correspondingly the degeneracy of 
the excited states is 1 as well). Now we use the fact that we 
have already proven our argument for the total number of Ma- 
joranas n’ = 8 — k, where the ground state degeneracy is 2, 
and the two ground states are different by a fermion particle 
number. 

Let us now bring the n’ Majoranas and 8m + k original 
ones together. The ground state of the decoupled system is 
one with an empty and an occupied fermion level. To cou- 
ple the two ground states one needs a degenerate empty and 
an occupied fermion level in the rest of the system, which 
is absent due to the ground state degeneracy 1. To prove 
this statement, one needs to rewrite the total Hamiltonian in 
the block-diagonal form, where one block corresponds to the 
empty fermion level, and another to the occupied. Then one 
notices that as each of the block is non-degenerate and as there 
is asymmetry of the total Hamiltonian under changing the oc- 
cupation of the fermion, the doubly-degenerate ground state 
stays doubly-degenerate even after coupling to the rest of the 
system. Therefore for 8m + k + n’ = 8(m + 1) Majoranas 
there is a region of parameter space, where the ground state 
is doubly degenerate, which contradicts the assumptions, and 
consequently proves the statement for k = 2, 6. 

We can now add to 8m + 2 Majoranas another two. We no- 
tice that there are two fermionic modes, empty and occupied, 
forming the ground state manifold without the coupling be- 
tween the two systems. It is clear that since there is no direct 
tunneling between the Majoranas, the only coupling allowed 
in the model consists of the two new Majoranas and the occu- 
pation number of the old fermion. Therefore we are left with 
the ground state degeneracy 2, but this degeneracy is between 
two states different by a bosonic operator, which proves the 
constraint on the ground state degeneracy for k = 4. 

Now we notice that k = 3,5 is different from k = 4 by a 
single Majorana fermion. By adding or subtracting the Ma- 
jorana it is impossible to remove the bosonic degeneracy of 
the ground state and the degeneracy should be enhanced by a 
factor \/2 over k = 4. This concludes the proof. 


IV. LAPA-TEO-HUGHES MODEL 


We now turn to the 1d models of interacting MZMs. Lapa, 
Teo and Hughes#4! (LTH) introduced a model of fermions in 
class BDI that is forced to be in the trivial phase by additional 
inversion symmetry in the absence of interactions. When in- 
teractions that respect the symmetries are introduced, the sys- 
tem, however, might transition to the topologically non-trivial 
phase characterized by gapped bulk and 4-fold ground state 
degeneracy associated with a pair of Kramers doublets bound 
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FIG. 3: (Color online) Lattice geometries for the interacting Majo- 
rana models. Color blue(red) indicates (anti)vortex Majorana modes, 
while solid and dashed lines indicate interaction and hopping terms, 
respectively. Panels a) and c) represent the 2-leg and the 4-leg LTH 
ladder, respectively. These two models are invariant under particle- 
hole (2), time reversal O, and inversion P operations. b) The 2-leg 
LTH model under the Majorana basis transformation, as discussed in 
the text. 


to the two edges of the 1D system. An experimental realiza- 
tion of the 1D LTH model, using MZMs in vortices and an- 
tivortices in the surface of a TI, has been proposed=21 The pro- 
posed setup is depicted in Fig. Bc) and consists of alternating 
clusters of 8 vortices and antivortices (that support 4-fermion 
interactions) connected to one another via tunneling terms. In 
the limit of strong interactions the ground state can be thought 
of as a direct product of unique ground states associated with 
each 8-vortex cluster. The ground state degeneracy is associ- 
ated with the quartet of MZMs located at each end of the chain 
and is protected by ©. This is a genuine interaction-enabled 
topological phase as no nontrivial phase can exist in a system 
with these symmetries in the absence of interactions“! Im- 
portantly, it persists as a stable phase for a finite strength of 
hopping t. 


In this Section we study the LTH model in detail by numer- 
ical techniques. We determine its phase diagram, ground state 
degeneracy as well as the entanglement spectrum. To facili- 
tate this study we introduce a closely related 1D model that is 
exactly soluble. This is depicted in Fig. Bia) and consists of 
alternating clusters of 4 vortices and antivortices. We call it 
the “2-leg LTH model”. The original “4-leg LTH” model can 
thus be thought of as two coupled 2-leg models. We also dis- 
cuss potential realizations of the LTH-type models in 2D and 
3D and their relevance to experimental systems. 


A. 2-leg LTH ladder 


The 2-leg LTH model Fig. Bia) can be regarded as com- 
posed of two parallel Majorana chains, which preserve time 
reversal © and the reflection symmetry. Since the presence 
of inversion symmetry trivializes the Fidkowski-Kitaev Zg 2- 
channel chain) this model should exhibit only topologically 
trivial phases even in the presence of interactions. Using the 
a, @ notation as illustrated in the right half of Fig. Bia) the 
Hamiltonian is written as 


M 2 
H = — it > S > (0ok-1,aF2k,a + Bb¢,a02k-+1,a) 
k=—M a=1 
M M-1 
+9 S> ho(Box, Box) +9 S> ho (Orn +1, K-41): 
k=—M k=—M 


(26) 


We consider a system with N’ Majorana sites along the chain, 
N’ being a multiple of 4. For convenience we define (half) in- 
teger valued MM = (N’ — 4)/8. Index & is then also (half) 
integer and extends between —M and M. The interaction 
ha(An: Xn) = AntXn,1An,2An,g: The hopping between the 
same type of Majoranas, which breaks the time reversal sym- 
metry, is absent in the Hamiltonian. The symmetries allow for 
interaction terms containing two a and two £ operators but we 
expect these to be small compared to the direct hopping terms 
already included in H. We thus include only 4-fermion terms 
between the same type of MZMs. 

The model is integrable because it has an extensive num- 
ber of constants of motion. Specifically, it is easy to see that 
products 


1 

AM = 091,10 9k—1,292k, 1 B2k,2 (27) 
2 

A?) = Bbx 1 Bb 200K41,1026-41,2 


commute with H for all k. In the following we introduce a 
non-local transformation that allows one to replace these op- 
erators by c-numbers and maps the problem onto a collection 
of non-interacting 1D Kitaev chains. 

To effect the transformation it is expedient to relabel the 
MZM operators as indicated in the left half of Fig. Bia). The 
Hamiltonian can be written in an economical way 


N 
| it ener + Yo1417214+2) 
1=0 
N-1 
-g S- Ya142V21437214212143> (28) 
1=0 


where NV = (N’ — 2)/2. 

The operators transform 7; > (—1)!y¥5n 43 j and y; > 
(—1)/72n+43~-; under the inversion operation, of which the 
inversion center is located at the middle of the (NV + 1)-th 
and (N + 2)-th sites since (a, a’, 3, 8’) + (a’, —a, 8’, —8). 
(The minus signs are coming from the spin-5 of the Fu-Kane 
model, see [IV D)) Since inversion symmetry in a 1D system 


plays an effective role of reflection symmetry, it will be shown 
in that the Hamiltonian is invariant under all of the sym- 
metries in class BDI+R__ S471 

Any 1D chain in this reflection symmetry class without in- 
teractions is always in the trivial phase. For this 2-leg model 
even in the presence of the interactions the system is still triv- 
ial, but it is exactly solvable. 

The transformation, analogous to the Jordan-Wigner trans- 
formation is 


k<j 
V5 =P, ( Il ia ) ssa, (29) 
k odd 
k<j 
Yj —— ip,( II ia) ibj-r4, (30) 
k even 
(31) 
where 
7 1, gj odd, 
si { —i, j even. (62) 


Here a; and b, are Majorana operators satisfying 


{ax, Qn} = 25k,k, {dk, bee} = 26k,Kr, Lak; de } = 0, 
(33) 


and /; runs from 0 to NV. It is straightforward to check that this 
transformation does not change the original +, +’ commuta- 
tion relations. Now we can rewrite the original Hamiltonian 
in the new Majorana basis. For odd 7 it holds 


VWi+1 =j-145 41, a 
Vi Vj-41 =tbj-145bj 4541, 2) 
while for even 7 
ViVi = — tbj-14;bj 0541, Ce) 
VGVj41 = — bj-145 41. G7) 


Therefore the Hamiltonian is rewritten as 


N 
H=— itS- boa2i42(1 + ta214162141) 
1=0 
N-1 
+9 > iar42ba142. (38) 
1=0 


We note that the number of Majoranas in a, b basis is two 
more than the number of Majoranas in +, 7’ basis as shown in 
Fig. Bib). The reason is that we have introduced two extra free 
Majoranas ao and bzn+2. These are necessary to satisfy the 
canonical commutation relations and the redundancy should 
be kept in mind when discussing the ground state degeneracy 
in the transformed basis. 

After the transformation, it is clear the system has many 


integrals of motion s; = 7a2;41b2;41. These are the same 
quantities as defined in Eq. (27). We can replace them by 
their respective eigenvalues s; = +1. In each sector, labeled 


by the set of quantum numbers {s;}, the transformed Hamilto- 
nian is non-interacting and represents a “broken” Kitaev chain 
illustrated in Fig. 3b). It has alternating bonds with hop- 
ping strength given by g and (1 + s;)t. In the sector with all 
8, = +1 this model describes a single Kitaev chain with a 
phase transition between its two gapped phases at t = g/2. 
In the sector with all s; = —1 the system consists of discon- 
nected monomers and the spectrum is a pair of flat bands; ¢ 
does not enter. In the mixed sector, where some hoppings 
are present and some are absent, we have an array of the dis- 
connected Kitaev chains of finite length. It is clear that the 
wider the bandwidth of the Hamiltonian, the lower the energy 
of the ground state in a given sector. Therefore, we expect 
the absolute ground state to be in the s; = +1 sector. We 
have confirmed this by explicit numerical simulation of the 
Hamiltonian (38) on a lattice up to 19 integrals of motion (20 
Majoranas in the hopping model, 58 in total). 

We now study the ground state s; = +1 sector. It is 
known that the Kitaev chain is dual to the transverse-field 
Ising model 48! The phase transition between its two gapped 
phases is therefore in the transverse-field Ising model univer- 
sality class. Indeed it is possible to transform the 2-leg LTH 
Hamiltonian onto a set of broken Ising chains, either directly 
going from y operators to spin-4 operators, or in two steps via 
the broken Kitaev chain Eq. (38). This also helps to under- 
stand the nature of the gapped phases of the original model. 
In the transformed basis for t > g/2 the absence of ground 
state degeneracy leads to the topologically trivial phase. For 
t < g/2 Majorana zero modes appear at each end in the a, b 
basis. When we assume N > 1, the form of the end modes is 
given by 


1 ya 
B= b 
oo (39) 
N l 
1 2t 
A= ihe 4 
a =) a2N—21 (40) 


The fermion operator B + Ai switches one ground state to the 
other. Although these Majorana operators are localized in a, b 
basis, the fact that the transformation is non-local means 
that in the original 7, 7’ basis these operators do not represent 
edge degrees of freedom. In the original basis t < g/2 is a 
conventional broken symmetry phase, akin to the ferromag- 
netic phase in the Ising model. [The symmetry that is sponta- 
neously broken in Hamiltonian for t < g/2 is generated 
by yj, —> —7;,- It remains unbroken in the other phase.] The 
2-leg LTH model undergoes a phase transition but both phases 
are topologically trivial as expected on the basis of general ar- 
guments presented above. 

To confirm the above conclusions we have performed 
DMRG calculations to compute the energy spectrum, the en- 
tanglement spectrum, and the central charge of the system as a 
function of dimensionless coupling g/t. The spectra in panels 
(a), (b) of Fig. [4|show that the phase transition indeed occurs 
at g/t = 2 as the gaps close. The ground state is doubly 
degenerate for g > 2¢ and the entanglement spectrum also 
shows 2-fold degeneracy in the symmetry broken phase. Im- 
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FIG. 4: (Color online) The numerical results of the 2-leg LTH model are computed using DMRG as g/t varies. Here N’ indicates the length 
of the chain, N’ = 2N + 2 in Eq. (28). a) AE; = E; — Ep is the energy difference between the ground state and the i-th excited state in a 
system with open boundary conditions. b) Entanglement spectrum for periodic boundary conditions. c) The central charge c. At the critical 
point it is well defined and size-independent. Furthermore, the energy gap closes and the number of the degenerate states in the entanglement 
spectrum changes at g/t = 2. Thus, the critical point is at g/t = 2 as expected based on the analytic solution. In panels (a) and (b), N’ = 120. 


portantly, the ground state degeneracy in this case occurs for 
both open and periodic boundary conditions, confirming that 
g > 2t is a conventional broken symmetry phase. To con- 
firm the phase transition property, we numerically compute 
the central charge c using. its relation with the entanglement 
entropy of ground states! 


c N' . 
- log(— sin 
3 T 


S(n) a7) + So. (41) 
Here N’ is the system size, n is the size of subsystem, S(n) 
is the corresponding entanglement entropy and So is a con- 
stant. At the critical point Fig. l4fc) indicates the value 4 of 
the central charge in agreement with our expectation that the 
phase transition is in the universality class of the transverse- 


field Ising model. 


B. 4-leg LTH ladder 


We proceed with the examination of the 4-leg LTH model 
whose geometry is illustrated in Fig. Bic). In the absence 
of interactions the restriction of inversion symmetry prohibits 
any topologically non-trivial state. Let us start with a non- 
interacting 4-leg setup to illustrate this. Due to the restric- 
tion of the symmetries, only hopping terms between a and 
{8 MZMs are present in the Hamiltonian. If we focus on the 
nearest neighbor such hoppings the Hamiltonian reads 


M 4 


Hnop = —t S- S > t(ab,-1,5Fak,j + Pak,j @2k-+1,9) (42) 


k=—M j=1 


where is a (half) integer. The system is formed by a set 
of Majorana dimers and, importantly, all the Majoranas in the 
system are hybridized by the hopping. Therefore, there are 
no edge states and the system is topologically trivial. To gen- 
erate the topological phase, we introduce reflection symme- 
try preserving interactions. The dominant interaction involves 


groups of four Majoranas of the same type, 


Ay (Ar) =9'A1,1A1,241,3A1,45 (43) 
ho (Ay, Az) =91 (Anat Ar 27,2 + An,3A1,31,41,4) 
+ 921 2A12A1,3A1,3 (44) 


where ;,; stand for either a or 8 MZM. As shown in Fig. Bic) 
the Hamiltonian of the open chain including the interactions 
can be written as 


Hiri =Hhop oF Hint, (45) 
where 
M 
Hint = S- [Fy (O,-1) + | (Box) + hy (Box) + Py (Orr, 41) 
k=—M 
M M-1 
+ So ho(Box, B56) + S2 ho(org+1s K-41): 
k=—M k=—M 
(46) 


To implement periodic boundary conditions, we add an ex- 
tra interaction term ha (jy 41, O41) to the Hamiltonian. 
The Hamiltonian is invariant under inversion symmetry oper- 
ation {a,j 5} > {1p 55, O-1,5—y} and (81,5, 8) 3 > 
{Boi 5—53 —B_1,5—;} since the Fu-Kane model is an intrinsic 
spin-5 system ay SL 

Let us first consider the extreme case when the hopping is 
off, t = 0. The Hamiltonian then consists of decoupled clus- 
ters each containing 8 (anti)vortex Majoranas. We may com- 
pute the ground state of each such cluster described by 


Hsub = hy (A) 1 hy (X’) =F ho (A, Nn) (47) 


Assume gi, g2, g’ are positive and define 4 complex fermions 
dja = (Ajq +iXj,q)/2 in each cluster. The many-body wave- 
function can be expressed in the fermion basis |njn2n374), 
where nq is the eigenvalue of the fermionic number opera- 
tor (Aja = Oh alga The unique ground state is then given 
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FIG. 5: (Color online) Numerical results for the 4-leg LTH model. We define the interaction strength g = gi = go = g’ and perform DMRG 
with the system size N’ = 8M + 4 = 120. a) AE; = FE; — Ep is the energy difference between the ground state and the i-th excited state 
for the 4-leg LTH with open boundary conditions. b) The entanglement spectrum of the ground state for the 4-leg LTH with periodic boundary 
conditions as a function of g/t. The LTH model exhibits 4-fold degeneracy for the ground state and the entanglement spectrum when g/t is 
larger that the phase transition point g./t ~ 0.83. c) The central charge as a function of g/t. 


by (|0000) — |1111))/./2 with energy —29, — go — g’. As 
a result, with the periodic boundary conditions the 4-leg LTH 
model has a unique ground state in the limit when t = 0. Since 
there is a gap to the lowest excited state we expect this non- 
degenerate ground state to survive for some range of non-zero 
Li: 

For open boundary conditions that respect the reflection 
symmetry, four a's M-1,a at the left end and four a2)441,4 at 
the right end have only interactions described by hj (a’_»,7_ 1) 
and h)(a2,;4,) respectively. It can be easily seen that four 
many-body Majorana operators constructed from the edge 
MZMs 


/ / / / 
O_omM—1,1%-2M-1,2) %-2M—-1,30-2M-1,4) 


9M 4+1,10%2M41,2> %2M41,302M4+1,4 


commute with the full Hamiltonian. This implies 4-fold de- 
generacy of the ground state associated with the ends of the 
chain. As argued in Ref.|44|the emergent time-reversal © acts 
anomalously in this degenerate subspace (such that ©? = —1) 
and the states can thus be viewed as two Kramers doublets. 
We therefore expect the degeneracy to be robust against any 
perturbations that do not break © and as long as the bulk re- 
mains gapped. 

We now study the general case with hopping ¢ turned on 
by means of DMRG. For simplicity and concreteness we take 
g = 9. = 92 = g’ but this is by no means essential. Fig. 
fa) shows 4 lowest energy states obtained by DMRG as a 
function of g/t for open boundary conditions. This indicates 
a phase transition at g./t ~ 0.83 to the state with a 4-fold 
degenerate ground state, in agreement with the above analy- 
sis. Importantly, no such degeneracy is observed for periodic 
boundary conditions which confirms the topological character 
of the interacting phase. 

An alternative way to confirm the non-triviality of the phase 
is to consider the degeneracy of the entanglement spectrum. 
To compute the entanglement spectrum, the LTH model with 
periodic boundary condition is separated by the two vertical 


cuts between horizontally closest (anti)vortices. In agreement 
with the ground state degeneracy of the open chain, the en- 
tanglement spectrum Fig. [5{b) exhibits four-fold degeneracy 
since each cut contributes two-fold degeneracy from the wave- 
function (|0000) — |1111)) /\/2 at the cut. 

Fig. [5{c) displays the central charge computed using 
DMRG as a function of g/t. Similar to the case of the 2-leg 
LTH model it saturates at a size-independent value of c = $ 
when g = ge. This suggests that the phase transition in the 
4-leg LTH model is in the transverse-field Ising universality 
class. Understanding more fully the microscopic origin of this 
transition is an interesting problem which we leave for future 
study. 

For topological crystalline insulators and superconductors, 
the bulk topology is commonly determined by the presence or 
absence of protected gapless modes at the boundary, that are 
invariant under a spatial symmetry. However, this is not the 
case for this interaction-enabled topological phase. Each end 
of the LTH model, which maps onto the other under inver- 
sion, is not invariant under inversion so the inversion sym- 
metry does not protect the Majorana end-modes. Instead, 
the LTH model inherits its topology from the class BDI in 
1D, which hosts stable Majorana boundary modes protected 
by time reversal symmetry and particle-hole symmetry. The 
presence of these Majorana boundary modes can be induced 
by symmetry-preserving interactions even if the phase would 
otherwise be trivial in the non-interacting system with the 
same symmetries. 


C. LTH-type models in dimensions 2 and 3 


In two dimensions the interaction-enabled phases are possi- 
ble in symmetry class BDI with an additional C rotation sym- 
metry. In the non-interacting case with weak invariants v, and 
Vy in x and y directions respectively from 1D class BDI, the 
symmetry requires vy, = vy and vy = —v,. The only solution 
iS Vz = Vy = O and all such systems must be topologically 
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FIG. 6: (Color online) Proposed minimal geometry for an LTH-type 
model in 2D. Red (blue) circles indicate MZMs of type a (3). When 
the interaction dominates, the Majorana system in the bulk is gapped. 
The edge states are either gapless or spontaneously symmetry bro- 
ken. 


trivial. In the interacting case, as before, the 1D classifica- 
tion changes to Zg and v, = vy, = 4 becomes another pos- 
sible solution of the above criterion indicating an interaction- 
enabled topological phase. A similar argument can be made 
for a 3D system with C4 rotation symmetry about three or- 
thogonal axes. 


The 2D model we propose is depicted in Fig. [6] Similar to 
the 1D case it consists of alternating clusters of 8 interacting 
MZMs connected to one another by hoppings. To understand 
the phase diagram of the model we notice that in the absence 
of interactions all the MZMs have counterparts connected by 
hoppings ¢, and therefore the ground state is unique with a gap 
2¢ to all excitations. In particular, the edge is also gapped. On 
the other hand, in the absence of hoppings and in the pres- 
ence of generic interactions within each cluster containing 8 
Majoranas the system is gapped in the bulk, as all the clusters 
have non-degenerate gapped ground state. At the same time 
we notice that the system has gapless flat-band edge states as- 
sociated with the decoupled quartets of MZMs that exists at 
the boundary. Each cluster of 4 Majoranas on the edge has 
doubly-degenerate ground state. These two ground states are 
different by two Majorana operators, and therefore are related 
by emergent eat reversal symmetry © that acts anomalously 
in this subspace“ such that ©? = —1. Thus Kramers theo- 
rem applies and the edge state must either remain gapless or 
spontaneously break the time-reversal symmetry even away 
from the strongly interacting limit, as long as the bulk remains 
gapped. In the latter case it is doubly-degenerate. Indeed we 
notice that if the hoppings are turned on, the edge spectrum 
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FIG. 7: (Color online) Proposed minimal geometry for an LTH-type 
model in 3D. Red (blue) spheres indicate MZMs of type a (8). a) 
The unit containing 8 MZMs of both types. b) The full 3D structure. 
Notice that surfaces consist of purely a-type MZMs. 


becomes gapped with a doubly degenerate ground state. This 
can be seen as follows. The hoppings will perturbatively gen- 
erate additional interactions between the neighboring quartets 
of MZMs along the edge, such as that denoted in Fig. [6] as 
gz. Effectively, then, the edge is described by a 1D interacting 
model discussed in Ref. In the thermodynamic limit the 
model has been shown to exhibit a gapped doubly-degenerate 
ground state. 

A similar construction can be given for a 3D system. A unit 
cell, containing 8 MZMs of each type is displayed in Fig.[7{a). 
When such unit cells are stacked to form a cubic lattice a 3D 
version of LTH-type model emerges, Fig. [7{b). The discus- 
sion of its phases parallels the discussion of the 2D case and 
we will not repeat it here. On that same basis we expect the 
model to exhibit two phases, one topologically trivial for weak 
interactions and one topological when interactions are strong. 
In the latter case the bulk once again is non-degenerate and 
gapped while the surfaces perpendicular to the Cartesian axes 
are anomalous in that they exhibit either gapless states pro- 
tected by © or spontaneously break © and are then doubly 
degenerate and gapped. 


D. Proposed experimental realizations of LTH-type models 


We now discuss potential realizations of the interaction- 
enabled phase for LTH-type models in dimensions | and 2. 
Our starting point is a thin film (or a flake) of an STI whose 
surfaces have been made superconducting by coating with a 
thin layer of an ordinary s-wave SC, such as Al or Nb. When 
a perpendicular magnetic field B is applied to this structure 
vortices are induced in the top surface while antivortices are 
induced in the bottom surface. MZMs bound to these are our 
basic ingredients. We now imagine that both surfaces are pat- 
terned with an array of holes in the SC layer as indicated in 
Fig. The holes serve two purposes: (i) they pin vortices 
in the desired positions and (ii) they remove the undesired 
low energy Caroli-Matricon-de Gennes states that would oth- 
erwise exist in the cores of vortices in the ordinary SC. This 
way, vortex cores will only exist in the surface state of the TI 
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FIG. 8: (Color online) Proposed physical realizations of the of the LTH-type models in 1d and 2D. Superconductivity is induced in two surfaces 
of a thin TI film. The SC coating is patterned to create an array of pinning sites for vortices which then form structures that approximate the 
model geometries shown in Figs. [6] and 3] a) A quasi-1D system respecting the inversion symmetry can be used to realize the 4-leg LTH 
model. b) A 2D system with C4 rotation symmetry that can be used to realized the LTH model in 2D. 


and will be located in the right spatial positions. When the 
chemical potential pz of the TI is tuned to the Dirac point and 
when the TI film is sufficiently thin then the hopping ampli- 
tude t will be appreciable only between MZMs located in the 
adjacent vortices and antivortices, indicated by green lines in 
Fig. MZMs in the same surface are then connected only 
through 4-fermion interactions. If we furthermore take into 
account only the dominant interaction terms that occur be- 
tween the quartets of MZMs that are closest together then it 
is easy to see how devices in Fig. [6] approximate the LTH- 
type models in 1D (panel a) and 2D (panel b). We note that 
by layering the structure in Fig. [8{b) it is possible to envision 
creating also the 3D model but in that case it is not clear how 
the pattern of hole might be fabricated to pin the vortices in 
the interior layers. 

We now show that the LTH models in this heterostruc- 
ture belong to reflection symmetry class BDI+R__“Jand a 
pair of (anti)vortices preserves inversion symmetry, which is 
equivalent to reflection symmetry in 1D. Consider two vor- 
tices located at positions (40,0). The vortices we take into 
account by specifying a position-dependent phase of the pair- 
ing function 


A = Aget(9+t 9), (48) 


where 64 = tan~* 4. The pairing function is invariant 
under inversion (x,y) — (—a2,—y) while the 2D reflection 
symmetry in any direction is broken. With A(R) = A(—R), 
the Hamiltonian Hp in Eq. is invariant under inversion 
operation 


P—'Hpx(—p, —R)P = Hrx(p,R), (49) 
where P = it°a”. Hence, we are able to construct 1D LTH 
models with the distribution of vortices and antivortices obey- 
ing A(R) = A(—R) in Fig. [8] (a). To check the topological 
classification, let the reflection symmetry operator be Hermi- 
tian R = o* and R anticommute with © and =. Refs. 
show that the topology of such a 1D non-interacting chain in 
class BDI+__ is always trivial. 


We denote a pair of two Majorana zero modes located re- 
spectively at two (anti)vortices (+b, 0) as ay, a’_,(8s, 6" 4) 
in| They transform to a’_,, —ay(8"_,, —6») under reflec- 
tion symmetry since ina spin-5 system Pp? =-1. 

Similarly, the setup of 2D LTH model can be designed as 
shown in Fig. [8](b). The only difference is that each dot now 
represents two MZMs. This can be arranged by designing 
two pinning holes to be placed close to one another or else 
following more closely the octagonal pattern design displayed 
schematically in Fig. (6| In any case the holes are arranged 
so that the C’, rotation symmetry is preserved. Furthermore, 
the symmetry has to be checked microscopically. Consider 
four vortices located at (0,0) and (0, £b). The phase of the 
pairing function coming from the four vortices is given by 


es Aget(9e.0t9-4,0+00,4+90,—a) (50) 
where 02,,y, = tan~' 4“. It is easy to check that the pair- 
ing function is invariant under C’', rotation symmetry opera- 
tion (x,y) + (y, —2). In addition, the Fu-Kane Hamiltonian 
is invariant under C4 


Cy" Hex ((y, =2), (Py; —pz))C4 = Hgx (a, y); (nsDy))- 
(51) 


where Cy = wa (1+%o0,). Hence, the design indicated in Fig. 


[8{b) obeys the Cy symmetry and provides a possible realiza- 


tion the 2D LTH model. 


V. DISCUSSION AND CONCLUSION 


Majorana zero modes bound to vortices in the SC surface of 
a topological insulator present a unique opportunity to study 
the effect of strong interactions in a fermionic system. This 
is because the kinetic energy of such fermions may be made 
to vanish by tuning a single parameter, chemical potential py 
of the underlying TI. The interactions between MZMs, even if 
nominally not very strong, become the dominant energy scale 
in the problem. We have suggested a specific setup to experi- 
mentally probe the effect of MZM interactions on the system’s 


ground and excited states. We have discussed how to model 
OD and 1D phases in this setup, and have proposed extensions 
to higher dimensions. Below we touch on the immediate ex- 
perimental relevance of the setup discussed. 

The ingredients for the proposed setups are all in place. Su- 
perconductivity has been successfully induced in the surface 
state of the 3D TI by multiple groups2"“ The ability to tune 
the chemical potential to the Dirac point as we require, has 
also been demonstrated 5454 Finally, vortices have been im- 
aged in such device and spectroscopic evidence for MZMs 
in the cores of such vortices has been reported 58] We note that 
in the latter experiment the chemical potential was far away 
from the Dirac point so this specific system would not work 
for the purpose we envision in this paper. However, given the 
rapid pace of the experimental progress, we expect this last 
obstacle to be overcome in the near future. 

Once the above hurdle has been surmounted not much ad- 
ditional fabrication is required for our OD Corbino ring pro- 
posal discussed in Sec. III. We envision testing the effect of in- 
teractions and the resulting Fidkowski-Kitaev Zg periodicity 
before attempting to engineer the more complex 1D and 2D 
structures. The OD system should provide a definitive signa- 
ture of the presence of the MZM interactions by changing the 
periodicity of the occurrence of the zero-bias peak in the tun- 
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neling conductance. The energy scales of few degrees Kelvin 
we estimated show that the observation is possible with the 
currently available technology. 


An ambitious, longer-term goal will be to engineer the 1D 
and 2D LTH-type models perhaps based on the designs out- 
lined in Fig. [8] using the same technology. We note that be- 
sides creating the desired vortex patterns through artificial pin- 
ning, it is possible to obtain them naturally as a Josephson 
vortex lattice in a junction between two superconductors on 
top of a TIA Probing the electron degrees of freedom in such 
systems is feasible Use neling spectroscopy as has been 
demonstrated recently. 
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